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Abstract 

The 7r + S + , 7T + H°, K + p, K + n, and scattering lengths are calculated in mixed-action 

Lattice QCD with domain-wall valence quarks on the asqtad-improved coarse MILC configurations 
at four light-quark masses, and at two light-quark masses on the fine MILC configurations. Heavy 
Baryon Chiral Perturbation Theory with two and three flavors of light quarks is used to perform the 
chiral extrapolations. To the order we work in the three-flavor chiral expansion, the kaon-baryon 
processes that we investigate show no signs of convergence. Using the two-flavor chiral expansion 
for extrapolation, the pion-hyperon scattering lengths are found to be a„-+s+ = —0.197 ±0.017 fm, 
and a n + s o = —0.098 ± 0.017 fm, where the comprehensive error includes statistical and systematic 
uncertainties. 
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I. INTRODUCTION 



Lattice QCD calculations of meson-meson interactions have yielded predictions for physical 
scattering lengths at the few percent level [H, 0, [3j]. Several reasons underlie this striking 
accuracy. Firstly, at the level of the lattice calculation, Euclidean-space correlation functions 
involving pseudoscalar mesons have signal/noise ratios 1 that do not degrade, or only slowly 
degrade with time. Therefore, highly accurate fits of both single- and multi-meson properties 
are possible with currently available supercomputer resources. Recent calculations of multi- 
meson interactions relevant for the study of pion and kaon condensation have been performed 
with up to twelve mesons interacting on a lattice [U, 0, 0] with no appreciable degradation of 
signal/noise with time. Secondly, and perhaps more importantly, QCD correlation functions 
involving Goldstone bosons are subject to powerful chiral symmetry constraints. Since 
current lattice calculations are carried out at unphysical quark masses, these constraints 
play an essential role in extrapolating the lattice data to the physical quark masses, as 
well as to the infinite volume, and continuum limits. Chiral perturbation theory (%-PT) 
is the optimal method for implementing QCD constraints due to chiral symmetry, and in 
essence, provides an expansion of low-energy S-matrix elements in quark masses and powers 
of momentum Q. 

In contrast to the purely mesonic sector, recent studies of baryon-baryon interactions, 
the paradigmatic nuclear physics process, have demonstrated the fundamental difficulty 
faced in making predictions for baryons and their interactions jl, 0]. Unlike the mesons, 
correlation functions involving baryons suffer an exponential degradation of signal/noise at 
large times 2 and therefore pose a fundamentally different kind of challenge in extracting 



signal from data [TJJ. Furthermore, while baryon interactions are constrained by QCD 
symmetries like chiral symmetry, the constraints are not nearly as powerful as when there 
is at least one pion or kaon in the initial or final state. For instance, there is no expectation 
that the baryon-baryon scattering lengths vanish in the chiral limit as they do in the purely 
mesonic sector. In nucleon-nucleon scattering, the s-wave interactions are actually enhanced 
due to the close proximity of a non-trivial fixed point of the renormalization group, which 
drives the scattering lengths to infinity, thus rendering the effective field theory description 
of the interaction highly non-perturbative 12| . 

Given the contrast in difficulty between the purely mesonic and purely baryonic sectors 
described above, it is clearly of great interest to perform a lattice QCD investigation of the 
simplest scattering process involving at least one baryon: meson-baryon scattering. While 
pion-nucleon scattering is the best-studied process, both theoretically and experimentally, 
its determination on the lattice is computationally prohibitive since it involves annihilation 
diagrams. At present only a few limiting cases that involve these diagrams are being inves- 
tigated [13]. Combining the lowest-lying SU(3) meson and baryon octets, one can form five 
meson-baryon elastic scattering processes that do not involve annihilation diagrams. Three 
of these involve kaons and therefore are, in principle, amenable to an SU (3) heavy-baryon 
X-PT (HB^-PT) analysis [14] for extrapolation. The remaining two processes involve pi- 



1 Here the signal is the Monte Carlo estimate of the quantum correlation function evaluated on the lattice, 
while the noise represents the statistical fluctuations in the correlation function. 

2 A recent high-statistics study of baryon correlation functions on anisotropic clover lattices has found that 
the exponential decay with time of signal/noise occurs only asymptotically in time, and therefore, the 



signal/noise problem in baryon correlation functions is not nearly as severe as previously thought 10 ]. 
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ons interacting with hyperons and therefore can be analyzed in conjunction with the kaon 
processes in SU (3) HB^-PT, or independently using SU(2) HBx-PT. 

Meson-baryon scattering has been developed to several non-trivial orders in the 577(3) 
HBx-PT expansion in Refs. (lR Jl6| . extending earlier work on kaon-nucleon scattering in 
Ref. [13]. A very-recent paper [l8[ has reconsidered the SU(3) HB^-PT results using a 
different regularization scheme, and also derived results for pion-hyperon scattering in the 
SU(2) HB%-PT expansion. These works make clear that the paucity of experimental data 
make it is very difficult to assess the convergence of the chiral expansion in the three-flavor 
case. Further, in the pion-hyperon system, the complete lack of experimental data precludes 
a separate analysis in the chiral two-flavor expansion. A lattice calculation of meson-baryon 
scattering analyzed using x-PT is therefore useful not only in making predictions for low- 
energy scattering at the physical point, but also for assessing the convergence of the chiral 
expansion for a range of quark masses at which present-day lattice calculations are being 
performed. 

Meson-baryon scattering is also of interest for several indirect reasons. The K~n interac- 
tion is important for the description of kaon condensation in the interior of neutron stars [19] , 
and meson-baryon interactions are essential input in determining the final-state interactions 
of various decays that are interesting for standard-model phenomenology (See Ref. [loj for 
an example). Finally, in determining baryon excited states on the lattice, it is clear that the 
energy levels that represent meson-baryon scattering on the finite-volume lattice must be 
resolved before progress can be made regarding the extraction of single-particle excitations. 

The experimental input to existing y-PT analyses of meson-baryon scattering is exten- 



sively discussed in Refs. [15|, [16|, llTj, [18|. Threshold pion-nucleon scattering information is 



taken from experiments with pionic hydrogen and deuterium |21|, |22|, and the kaon-nucleon 



scattering lengths are taken from model-dependent extractions from kaon-nucleon scattering 



data [23j. There is essentially no experimental information available on the pion-hyperon 
and kaon-hyperon scattering lengths. There have been two quenched lattice QCD studies 
of meson-baryon scattering parameters: the pioneering work of Ref. [IiJ calculated pion- 
nucleon and kaon-nucleon scattering lengths at heayy pion masses without any serious at- 



tempt to extrapolate to the physical point, and Ref. [25| calculated the 7 = 1 KN scattering 
length and found a result consistent with the current algebra prediction. 

In this work we calculate the lowest-lying energy levels for five meson-baryon processes 
that have no annihilation diagrams: 7r + S + , 7r + S°, K + p, K + n, and 7\"°S° in a mixed-action 
Lattice QCD calculation with domain-wall valence quarks on the asqtad-improved coarse 
MILC configurations with b ~ 0.125 fm at four light-quark masses (m, ~ 291, 352, 491 
and 591 MeV), and at two light quark masses (m^ ~ 320 and 441 MeV) on the fine MILC 
configurations with b ~ 0.09 fm, with substantially less statistics on the fine ensembles. 
We extract the s-wave scattering lengths from the two-particle energies, and analyze the 
five processes using 577(3) HBx-PT. We find a rather conclusive lack of convergence in the 
three-flavor chiral expansion. We then consider 7r + E + and 7r + S° using SU(2) HB%-PT and 
find that we are able to make reliable predictions of the scattering lengths at the physical 
point. We find 

a w+s+ = -0.197 ±0.017 fm ; (1) 
a^+so = -0.098 ± 0.017 fm , (2) 

where the errors encompass statistical and systematic uncertainties. The leading order x _ PT 
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(current algebra) predictions for the scattering lengths are given by [26]: 

a w+s+ = -0.2294 fm ; (3) 
a^+so = -0.1158 fm . (4) 

Ultimately, either the chiral extrapolation should be performed after a continuum limit 
has been taken, or one should use the mixed-action extension of HB%-PT to perform the 



chiral extrapolations [271 . |28| . However, our results on the fine MILC configurations are 
statistics-limited and not yet sufficiently accurate to make this a useful exercise. Further, 
the explicit extrapolation formulas for the meson-baryon scattering lengths have not yet 
been determined in mixed-action X"PT. Despite these limitations, we expect the corrections 
from finite lattice spacing to be small for two principle reasons. Firstly, the meson-baryon 
scattering lengths are protected by chiral symmetry and therefore the (approximate) chiral 
symmetry of the domain wall valence fermions used in this work protects the scattering 
lengths from additive renormalization, which can be explicitly seen in the construction of 



the mixed-action baryon Lagrangian in Ref . |28J . The mixed-action corrections do not appear 



until next-to-next-to leading order in the chiral expansion of the meson-baryon scattering 
lengths. Secondly, our previous experience with this mixed-action lattice QCD program leads 
us to expect that discretization effects will be well-encompassed within the overall errors we 
quote. In our precise calculation of meson-meson scattering, the predicted mixed-action 
corrections 0, 30] were smaller than the uncertainties on a given ensemble [l|, Hj|. 



This paper is organized as follows. In section [TT] we isolate the five meson-baryon pro- 
cesses with no annihilation diagrams that are calculated in this work. We briefly review the 
standard Luscher method for extracting the scattering amplitude from two-particle energy 
levels in a finite volume in section IIII1 Particulars regarding the mixed-action lattice cal- 
culation and fitting methods are provided in section [TV] Additional details can be found 



in Ref. [311 ] . Mixing between two of the meson-baryon channels with the same quantum 
numbers is discussed in section |V] In section I VI I we consider chiral extrapolations of the 
lattice data using SU(3) HBx-PT, and in section IVHI we analyze the pion-hyperon lattice 
data using SU(2) HBx-PT. Finally, we conclude in section IVIIII 



II. MESON-BARYON SCATTERING PROCESSES 

It is a straightforward exercise to construct the six scattering channels involving the lowest- 
lying octet mesons and baryons that do not have annihilation diagrams, and to determine 
their isospin. 3 The particle content, isospin, and valence quark content of these meson- 
baryon states are shown in Table [Q We adopt the notation of Ref. [15|], denoting the 
threshold T-matrix in the isospin basis as T^, where / is the isospin of the meson-baryon 
combination, is the meson, and B is the baryon. The five elastic meson-baryon scattering 
processes that we consider are then in correspondence with the isospin amplitudes according 



3 The 7T + S° and K°Y, + systems have the same quantum numbers, and therefore require a mixed channel 
analysis in order to extract the K°Yj + scattering length. This is discussed in Section fVl 
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Particles 


Isospin 


Quark Content 


7r+x; + 


2 


uuuds 


-7T+- 

7T i i 


Of z 


tlttClbO 


K + p 


1 


uuuds 


K+n 


and 1 


uudds 


K°Z+ 


3/2 


uudss 


K°E° 


1 


udsss 



TABLE I: Particle content, isospin, and valence quark structure of the meson-baryon 
states calculated in this work. As is clear from the valence quark content, these meson- 
baryon states have no annihilation diagrams. 



to 



Tl rp(2) rp rp{3/2) 



rp _ rp(l) . T , I'TW _l_ T^ ) \ ■ T 1 — 

l-K+p — J-KNi J-K+n — \ KN KN/ ' i K S° — 



These threshold T-matrices are related to the scattering lengths a^s through 



(5) 



T^b = 4vr [l + ^ J a^B , (6) 
where m$ is the meson mass and mg is the baryon mass. 

III. FINITE- VOLUME CALCULATION OF SCATTERING AMPLITUDES 

The s-wave scattering amplitude for two particles below inelastic thresholds can be deter- 
mined using Liischer's method [13], which entails a measurement of one or more energy 
levels of the two-particle system in a finite volume. For two particles with masses and 
tub in an s-wave, with zero total three momentum, and in a finite volume, the difference 
between the energy levels and those of two non-interacting particles can be related to the 
inverse scattering amplitude via the eigenvalue equation [32J 

pcoa(p) = \s(^) , (7) 



ttL \2n / 

where 8{p) is the elastic-scattering phase shift, and the regulated three-dimensional sum is 

|j|<A 



S (^) = - 47rA • ( g ) 



"j \j\ -v 



The sum in Eq. (JSJ) is over all triplets of integers j such that |j| < A and the limit A — > oo 
is implicit [33| . This definition is equivalent to the analytic continuation of zeta-f unctions 
presented by Liischer 32|. In Eq. (j7]), L is the length of the spatial dimension in a cubically- 



symmetric lattice. The energy eigenvalue, E n , and its deviation from the sum of the rest 
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masses of the particle, AE n , are related to the center-of-mass momentum p n , a solution of 
Eq. ©, by 



AE n = E n - - m B 
+ ... , 



Pn + m 6 + V Pn + m B ~ m <t> ~ m B 



Pi 



(9) 



where fi^g is the reduced mass of the meson-baryon system. In the absence of interactions 
between the particles, |p cot <5| = oo, and the energy levels occur at momenta p = 27rj/L, 
corresponding to single-particle modes in a cubic cavity with periodic boundary conditions. 
Expanding Eq. (J7J) about zero momenta, p ~ 0, one obtains the familiar relation 4 



AEn 



2na 



^bL 3 



1 + ci- + c 2 - 



with 



|j|<A 



Cl 



- V ^7 - 4A = -2.837297 
and a is the scattering length, defined by 



c 2 



(10) 



6.375183 ,(11) 



lim 



tan 5 (p) 
P 



(12) 



As the finite-volume lattice calculation cannot achieve p = (except in the absence of 
interactions), in quoting a lattice value for the scattering length extracted from the ground- 
state energy level, it is important to determine the error associated with higher-order range 
corrections. 



IV. LATTICE CALCULATION AND DATA ANALYSIS 

In calculating the meson-baryon scattering lengths, the mixed-action lattice QCD scheme 



was used in which domain-wall quark (33.1351. 36 . 37. |38| propagators are generated from a 



smeared source on n/ = 2 + 1 asqtad- improved |39l. |40|| rooted, staggered sea quarks |4l| 
To improve the chiral symmetry properties of the domain- wall quarks, hypercubic-smearing 
(HYP-smearing) 



43l . 1441 ] was used in the gauge links of the valence-quark action. In 
the sea-quark sector, there has been significant debate regarding the validity of taking the 

"IE EE HE 



fourth root of the staggered fermion determinant at finite lattice spacing [4 



50, £11, 52 



54, 55, 56, 57, 58 



While there is no proof, there are arguments to suggest 
that taking the fourth root of the fermion determinant recovers the contribution from a 
single Dirac fermion. The results of this paper assume that the fourth-root trick recovers 
the correct continuum limit of QCD. 



4 In order to be consistent with the meson-baryon literature, we have chosen to use the "particle physics" 
definition of the scattering length, as opposed to the "nuclear physics" definition, which is opposite in 
sign. 
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The present calculations were performed predominantly with the coarse MILC lattices 
with a lattice spacing of b ~ 0.125 fm, and a spatial extent of L ~ 2.5 fm. On these 
configurations, the strange quark was held fixed near its physical value while the degenerate 
light quarks were varied over a range of masses corresponding to the pion masses shown in 
Table [II] See Ref. jUJ for further details. Results were also obtained on a coarse MILC 
ensemble with a spatial extent of L ~ 3.5 fm. However, this data is statistics limited. In 
addition, calculations were performed on two fine MILC ensembles at L ~ 2.5 fm with 
b ~ 0.09 fm. On the coarse MILC lattices, Dirichlet boundary conditions were implemented 
to reduce the original time extent of 64 down to 32, which saved a nominal factor of two in 
computational time. While this procedure leads to minimal degradation of a nucleon signal, 
it does limit the number of time slices available for fitting meson properties. By contrast, 
on the fine MILC ensembles, anti-periodic boundary conditions were implemented and all 
time slices are available. 





Ensemble 


m n (MeV) 


bm\ 


bm s 


7 dw f 

bra^ 


i dw f 

bm s 


10 3 x bm res a 


# of props 


(i) 2064f21b676m007m050 


291 


0.007 


0.050 


0.0081 


0.081 


1.604 ±0.038 


1039 x 24 


(u) 


2064f21b676m010m050 


352 


0.010 


0.050 


0.0138 


0.081 


1.552 ±0.027 


769 x 24 


(Hi) 


i 2064f21b679m020m050 


491 


0.020 


0.050 


0.0313 


0.081 


1.239 ±0.028 


486 x 24 


(iv) 


2064f21b681m030m050 


591 


0.030 


0.050 


0.0478 


0.081 


0.982 ± 0.030 


564 x 24 


(v) 


2864f21b676m010m050 


352 


0.010 


0.050 


0.0138 


0.081 


1.552 ±0.027 


128 x 8 


(vi) 


2896f21b709m0062m031 


320 


0.0062 


0.031 


0.0080 


0.0423 


0.380 ± 0.006 


1001 x 8 


(vii) 


2896f21b709m0124m031 


441 


0.0124 


0.031 


0.0080 


0.0423 


0.380 ± 0.006 


513 x 3 



"Computed by the LHP collaboration for the coarse ensembles. 



TABLE II: The parameters of the MILC gauge configurations and domain-wall propaga- 
tors used in this work. The subscript I denotes light quark (up and down), and s denotes 
the strange quark. The superscript dwf denotes the bare-quark mass for the domain-wall 
fermion propagator calculation. The last column is the number of configurations times the 
number of sources per configuration. Ensembles (i)-(iv) have L ~ 2.5 fm and b ~ 0.125 fm; 
Ensemble (v) has L ~ 3.5 fm and b ~ 0.125 fm; Ensembles (vi),(vii) have L ~ 2.5 fm and 
b ~ 0.09 fm. 

The correlation function that projects onto the zero momentum state for the meson- 
baryon system is 

CM*) =^ ^(0 t (t,x)^(t,y)0(O, 0)5,(0,0)) , (13) 

where Vij is a positive-energy projector. For instance, in the case of K + p, the interpolating 
operators for the K + and the proton are 

0(t,x) = K + (t,x) = s(t,x)7 5 tt(t,x) ; 

Bi(t,x) = pi(t,x) = e a6c <(t,x) ( M fcT (t,x)C 75 rf c (t,x)) . (14) 

The masses of the mesons and baryons are extracted using the assumed form of the large- 
time behavior of the single particle correlators as a function of time. As t — > oo, the ground 
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state dominates; however, fluctuations of the correlator increase with respect to the ground 
state. The meson and baryon two-point correlators, C^(t) and C B (t), behave as 

C^t) - Ax e- m * \ C B (t) - A 2 e~ mB 1 , (15) 

respectively, in the limits t — > oo and L — > oo. In relatively large lattice volumes the 
energy difference between the interacting and non-interacting meson-baryon states is a small 
fraction of the total energy, which is dominated by the masses of the mesons and baryons [l[ . 
In order to extract this energy difference the ratio of correlation functions, G^ B (t), is formed 

GMt) - = X>» e ~ AEn 1 ' (16) 

where AE = AEq is the desired energy shift. With AE, and the extracted masses of the 
meson and baryon, the scattering length can be calculated using Eqs. (JTj) and (jHJ), or, if 
a << L, from Eq. ([TO]) . For the meson-baryon scattering lengths calculated in this work, 
the difference between the exact and perturbative eigen-equations is negligible. 

A variety of fitting methods have been used, including standard chi-square minimization 
fits to one and two exponentials. Generalized effective energy plots are particularly useful 
for analyzing the lattice data and for estimating systematic errors [Io|. These plots are 
constructed by taking the ratio of the correlators at times t, and t + nj (where n j is an 
integer) 

< B = -log L C **® , ) , AEf B = -llog ( G yW ) . (17) 

With rij — 1, the standard effective mass and energy plots are recovered. Generalized 
effective masses form a system of linear equations for each rij over the time interval where 
the data is fit. For instance, if the interval is given by At = t2—ti, then there is one equation 
for m eff at each t, for any rij that fits within At. The equations can be solved for m eff by 
casting them into the form of the so-called normal equation [59]. Since each rij constitutes 
a different effective mass plot, the number of degrees of freedom is increased significantly. 
This method provides a fitting routine that is faster than standard least-squares fitting. 
Additional details regarding the utility of generalized effective mass and energy plots can be 
found in Ref. [60l ] . 

The interpolating operator at the source is constructed from gauge-invariantly-smeared 
quark field operators, while at the sink, the interpolating operator is constructed from either 
local quark field operators, or from the same smeared quark field operators used at the 
source, leading to two sets of correlation functions. For brevity, we refer to the two sets 
of correlation functions that result from these source and sink operators as smeared-point 
(SP) and smeared- smeared (SS) correlation functions, respectively. By forming a linear 
combination of the SP and SS correlation functions, — aC^ SF \ we are able to remove 

the first excited state, thus gaining early time slices for fitting [6(|. This effect is illustrated 
in Fig. [H which is the effective AE n +-£+ plot for coarse MILC ensemble {ii). We plot C (ss) , 
C^ SP \ and — aC' SP ' with a tuned to remove the first excited state. The effective 

energies, effective masses, and energy splittings are plotted for coarse MILC ensemble (ii) 
in Figs. [21 [31 and [H All of the necessary quantities needed for extraction of the scattering 
lengths are contained in Table IHH which also contains the sum of meson and baryon masses 
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«Y C (SS, -aC <SP) 




FIG. 1: Effective AE 7T +y,+ plot for coarse MILC ensemble (ii) from correlation functions 
C( ss \ C*( SP ) and — aC^ SP \ By taking the linear combination with a tuned to re- 

move the first excited state, earlier time slices are gained for fitting. 



at each quark mass. Fig. [5] shows the results for all five processes, and the behavior of 
Eq. ([7]), versus the interaction energy, presented in terms of the dimensionless quantities 
p cot 5 1 m-x and AE / m n . The curve shown in Fig. is p cot 5 / m n for the case of = itik , 
and tub = Tn p , as AE/m^ is varied. S(r/) in Eq. (jSJ) is a function of the meson and baryon 
masses, so there will be a unique curve for each combination of and m#. Consequently, 
the K + p, and K + n data points fall on this curve. 



V. THE MIXED CHANNEL 

As is clear from Table I, the 7r + S° and K°Y, + states carry the same global quantum 
numbers, and therefore couple to the same energy-eigenstates in the finite lattice volume. For 
energies above both kinematic thresholds, a determination of the three scattering parameters 
associated with these states (two phases and one mixing-angle) requires a coupled-channel 
analysis. Therefore, three energy levels above both kinematic thresholds must be determined 
in the lattice calculation to fully characterize scattering in this kinematic regime. In the 
present lattice volumes, the two-particle energies in these channels are close to the respective 
kinematic thresholds, and the energy of the lower-lying 7r + H° state (which is below the A'°E + 
threshold) is determined by the low-energy elastic scattering parameters, making it amenable 
to analysis using Eqs. (|7|h (jgjh P and (fT0|). 

A priori, one would expect both the 7r + S° and K°T, + interpolating operators to couple to 
a common ground state (dominantly the 7r + S° state), with a ii" E + -related level as the first 
excited state (for the lattice volumes considered here, the non-interacting 7r + S° system with 
two units of relative momentum has an energy considerably above the i^°E + threshold). 
Interestingly, within our statistical and systematic uncertainties, we find distinct energy 
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7t + E + Total EmO 10 



7i + S° Total E mOlO 




10 12 14 16 



(a) 



1.16 



1.14 



P-l 



1.12 



1.1 



1.08 



m +m_=1.124 
n J = 2 



n 5 h = n 5 5 



t/b 



10 12 14 16 



(b) 



1.16 



1.14 



1.12 



1.1 



1.08 



K p Total E mOlO 



nv+m =1.11 1 



5 5^55 



II 



4 6 



t/b 



(c) 



10 12 14 



16 



K n Total EmO 10 




1.08 - 



(d) 



13 



1.25 



1.2 



1.15 



K°E + Total EmO 10 



m+m =1.231 
m +ml=1.124 



-^-n-^- n -^-5-- 5 --|-^-|- 



10 



t/b 



12 14 



16 



(e) 



K°S° Total EmO 10 



1.34 




(f) 



FIG. 2: Effective energy plots of the six meson-baryon processes shown in Table [H The 
plots are from MILC ensemble (ii), nj = 2, and the linear combination — aC' SP ' 

is plotted. The dashed line is the sum of the meson and baryon masses for each process, 
while the error bars represent the jackknife uncertainty Note that the bE axis of (e) is a 
factor of two larger in span than the other plots to encompass the dashed line at + 
m= = 1.124. 
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0.227 
0.226 
0.225 
0.224 
3 0.223 
0.222 
0.221 
0.22 



71 m010 



X per d.o.f= 0.683 
bm = 0.22305± 0.00025 
n=2 



U { I * * { 



10 



t/b 



(a) 



12 



14 



K mOlO 



a 



X per d.o.f= 0.349 
bm = 0.37816± 0.00026 
n =2 



5 5 



U j 5 * * j 



t/b 



10 12 14 



(b) 



0.76 



0.75 - 



0.73 



pmOlO 



X per d.o.f = 0.44 

bm = 0.73242± 0.00306 

n J = 2 



* * li [ \\\ 



8 10 12 

t/b 



(c) 



14 



16 



0.87 



E 



0.85 



0.84 



Z mOlO 



X per d.o.f = 1.592 
bm = 0.85305± 0.00191 



5 5 J 



(d) 



8 10 12 14 

t/b 



16 



E°m010 



0.92 



| 0-9 



0.89 



X'per d.o.f = 0.814 
bm = 0.900910.00134 
n J = 2 

E 5 H I 1 1 I I I 



10 12 14 16 



t/b 

(e) 

FIG. 3: Single particle effective mass plots for coarse MILC ensemble (it). Here we choose 
nj = 2, and the linear combination — aC^^ is plotted. The inner shaded bands are 
the jackknife uncertainties of the fits to the effective masses, and the outer bands are the 
jackknife uncertainty and systematic uncertainty added in quadrature over the indicated 
window of time slices. 
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FIG. 4: Meson-baryon effective energy difference plots for coarse MILC ensemble (ii). 
Here we choose n,j = 2, and the linear combination — aC' SP ' is plotted. The inner 
shaded bands are the jackknife uncertainties of the fits to the effective energy differences, 
and the outer bands are the jackknife uncertainty and systematic uncertainty added in 
quadrature over the indicated window of time slices. 
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Quantity 


m007 (z) 


mOlO (««) 


mOzO (mj 


m030 (zu) 


m T 


0.l8384(3l)(03) 


0.22305(25)(08) 


0.3l03l(38)(95) 


0.37513(44)(13) 


m k 


0.36783(32)(42) 


0.37816(26)(11) 


0.40510(33)(37) 


0.43091(66)(16) 


m p 


0.6978(6l)(08) 


0.7324(31) (10) 


0.8069(22)(14) 


0.8741(16)(05) 


m 2 


0.8390(22)(03) 


r\ oroi i 1 n\ i 'no\ 

0.8531(19)(08) 


0.8830(18)(17) 


n noi 0/1 \ /no\ 

0.9213(13)(03) 


ms 


n 00^0/1 o\/i p\ 

0.8872(13)(16) 


0.9009(13)(10) 


0.9233(18)(04) 


n r\ A r' ~i / 1 ,1 \ /no \ 

0.9461(14) (08) 


r 

U 


0.09257(16) 


0.09b00(14) 


0.10208(14) 


0.107b3(32) 


£ 

JK 


0.10734(10) 


0.10781(18) 


n 1 nn'7/ 1 / 1 t\ 

0.10976(17) 


n 1 1 oco/01 \ 

0.11253(31) 




0.0150(14)(08) 


0.0148(08)(13) 


0.0111(10)(08) 


0.0100(10)(11) 




0.00646(64) (98) 


0.00b2(05)(12) 


n n n i 1 ZccA/'.-ioA 

0.00431(b8)(43) 


0.00421(7b)(60) 




0.0140(22)(30) 


0.014b(15)(13) 


0.0092(10)(51) 


0.0087(16)(16) 


A EP 


0.0057(18)(lb) 


n nntr 1 /1 /I \ / nn\ 

0.0051(14)(09) 


0.003b(09)(12) 


0.0028(10)(11) 


A 77 


n n 1 1 o/no\/i \ 

0.0118(08)(13) 


n m or /nr \ / 1 /i \ 

0.0125(05) (14J 


n nnocr/no\/oi \ 

0.0085(08)(31) 


0.0086(16)(16) 




-2.12(16)(09) 


-2.36(09)(15) 


-2.30(15)(13) 


-2.36(18)(19) 


Q.7TS 


-1.08(09)(14) 


-1.19(09)(20) 


1 no/1 r \ / nn\ 

-1.08(15)(09) 


1 on / 1 o \ /1 c\ 

-1.20(18)(15) 




-2.8U(32)(44) 


-2.95(21)(19) 


/" n o \ /1 n\ 

-2.3(0.2)(1.0) 


ot/oi \ / \ 

-2.27(31)(32) 




-1.41(37)(34) 


-1.33(30)(21) 


-1.05(22)(30) 


n on/o'7\ / 01 \ 

-0.89(27)(31) 




-2.62(13)(21) 


-2.77(08)(23) 


10/1 r\// i o\ 

-2.18(15)(63) 


on/onN/ooN 

-2.29(30)(32) 


+ m p 


0.8817(61) 


n nrrr/oi \ 

0.9555(31) 


1.1172(23) 


1 1 no / 1 \ 

1.2492(18) 


7T7, + 777,v 


1.0229(23) 


1.0761(20) 


1.1933(19) 


1.2964(15) 


+ m= 


1.0710(14) 


1.1240(14) 


1.2336(19) 


1.3212(16) 


fflK + m p 


1.0657(61) 


1.1106(31) 


1.2119(23) 


1.3050(19) 


rriK + "is 


1.2069(23) 


1.2312(20) 


1.2881(19) 


1.3522(16) 




1.2550(14) 


1.2791(15) 


1.3284(19) 


1.3770(17) 



TABLE III: Lattice calculation results from the four coarse MILC ensembles which enter 
the analysis of the meson-baryon scattering lengths. The first uncertainty is statistical and 
the second uncertainty is systematic due to fitting. All quantities are in lattice units. 



levels from the two interpolating operators. This is consistent with strong coupling to the 
color-singlet constituents of the interpolating operator and only very weak couplings to 
states that require color rearrangement (see Fig. [2]). While this is suggestive that mixing 
between the states is small, a definitive interpretation requires an extraction of three energy 
levels above the kinematic thresholds of the tt + S° and and below the next kinematic 

threshold, in order to determine the three scatte ring parameters. The optimal way to extract 
these levels is to use the variational method 6l|, |62j, which requires the full matrix of 
correlation functions to be calculated, and diagonalized. The extraction of the scattering 
parameters would th en p roceed via an extension of the variational method to the coupled- 
channel scenario I 1 



64 



Due to our incomplete knowledge of the three mixed-channel energy levels, we do not 
attempt to extract any K°T, + scattering parameters in this work. 
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FIG. 5: pcotS/m^ versus AE^b/ 771 ^ f° r the five elastic scattering processes from coarse 
MILC ensemble (ii). The curve shown is pcotS/m n for the case of = mx, and mg = 



m p . 



VI. SU(3) HBxPT EXTRAPOLATION 

A. Scattering Length Formulas 

The scattering lengths of the five meson-baryon processes listed in Eq. (JSJ) are, to 0(m%) in 
SU(3) HBx-PT EQ, 



«7T+E + 



1 m s 
47r m,r + ms 



J 77 J 7V J 77 



(18) 



47r m^ + m~- 



J 77 J 77 J 77 



(19) 



47r m^: + 77i jv 



2m;c 2m 



/I 



Jk 



in 



A 



Pk 



(20) 



O'K+n 



m N 



47T nix + rriN 



rriK Tn 



Pk 



+ -#<?01 + ^jr+nM + 8/ii(/i) 
./A 



A" 



Pk 



(21) 



1 m= 

4-7T mx + m= 



2mA- 2m 



/ 2 



A' 



+ -7T LC 'i + ^7?o H o(/i) + 8/1123 (A*) 



rn 



A 



iA 



(22) 
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where we have defined Cqi = Cq + C\ and = /ii — /i2 + /13, and the loop functions are 
given by 

:w(aO = 2^1-7^^-2111 — -In — 

1 



2 9 % , T 

mi — mr. arccos 1 — 

rriK 2 



3F 2 m 7T D 2 m r 



(23) 



, n K J (3 m n m K 
3^+5° M = 4^4 1 -m 7r l--21n In 



/i 



— m 2 I 7r + arccos 



m„ 



7T 

+ 4 



3(£> - F) 2 m» — -(D + 3F) 2 m, 



m x 



(24) 



m 



7717, 



— + 31n^ 



in 



K 



ml In 



m K + a/ m|- - ml 



m„ 



mi 



m K arccos ■ 



m K 



m r 



71 



--(D-3F) 



2(D + F) 



m; 



mjj + m w 



+ (D + 5F)m ri 



(25) 



yK+n{v) 



m K 



m n 
In — 



In 



7T 



(D - 3F) 



(D + F) 



mz 



m K 
H 
1 



+ 



m 



A 



ml In 



m K + a/ 



A 



m„ 



+ ( 7jD + 3F ) m 

771,, + 771^ 



mi 



(26) 



A" 



A« 



TTiA 



, m w , ins- , m„ 
3 + 2 In — + In — + 3 In —2 



// 



+2\ m 



A" 



m£ m 1 3 



m„ 



m^- arccos • 



m K 



m r . 



(D + 3F) 



2(D-F)- 



mZ 



+ (D — 5F)m v 



(27) 



m v + m^ 

In what follows, we choose ji = A x = 4:7if n and evaluate f n at its lattice physical value j65[ , 
and we take m v from the Gell-Mann-Okubo formula. These choices modify the chiral ex- 
pansion at 0(m*) and are therefore consistent to the order we are working. The first 
mixed-action modification to these HB^-PT extrapolation formulas appear as corrections 
to these loop functions, y^B, and to the corresponding counterterms which absorb the scale 
dependence. Some of the mesons propagating in the loops appear as mixed valence-sea 
combinations, and thus the corresponding meson masses appearing in these functions are 
heavier by a known amount {66[. The precise form of the predicted corrections require a 
computation of the scattering processes with mixed-action/partially quenched x _ PT. 

Our physical parameters are consistent with Ref. [l8| (note that our decay constant 
convention differs by \/2). Namely, f n = 130.7 MeV, m n = 139.57 MeV, f K = 159.8 MeV, 
m K = 493.68 MeV, m N = 938 MeV, m s = 1192 MeV and m s = 1314 MeV. The axial 
couplings, D and F, for coarse MILC ensembles (ii)-(iv) are taken from the mixed-action 
calculation of Ref. [67|, and we extrapolate for coarse MILC ensemble («) using these values. 
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B. Extrapolation to the Physical Point 



For the purposes of fitting and visualization, it is useful to construct from the scattering 
lengths the functions r^ 1 ' 2 ) which are polynomials in m^. For the 7r + S + , K + p, and 
processes one defines 5 



,(i) _ 27ra/ f m<j} 



r^ = i£ M+^t =1; (28) 



r« M = -^S(l + ^) = l- Cim ,; (29) 



V m B / 2m 



123 (A x )m 2 , (30) 



and for the 7r + S°, and K + n processes one defines 



,( 2 ) __ ^afl ( ^ m^, 



r^ = ^1 + ^ =1; (31) 



r<?» --^(i + ^) = i-c„ im ,; (32) 



rgko = -— f 1 + — ) + — ^b(A x ) = 1 - C 01 m, - S/^AxK • ( 33 ) 

Notice that the left-hand sides of these equations are given entirely in terms of lattice- 
determined quantities, all evaluated under Jackknife, whereas the right-hand side provides 
a convenient polynomial fitting function. Plots of T^lo formed from the lattice data (all 
ensembles listed in Table HT|) versus the Goldstone masses are given in Fig. [61 We see evidence 
in this plot that the fine and large-volume coarse data are statistically limited as compared 
to the coarse data. Therefore, we include only the coarse data in our fits. The fine data is, 
however, indicative that lattice-spacing effects are small. 

In the three-flavor chiral expansion, we have an overdetermined system at both NLO 
and NNLO. While there are five observables, there are two Low Energy Constants (LECs) 
at NLO, Co and Coi, and two LECs at NNLO, h\ and h^- Fits of the LECs from each 
process at NLO are given in Table HVl and the corresponding values of the scattering lengths 
are given in Table [V] At NLO, the LECs are of natural size, and provide a consistent 
extraction within uncertainties. Correspondingly, the scattering lengths appear to deviate 
perturbatively from the LO values. The perturbative behavior of the scattering lengths at 
NLO is evident from the plots of T NLO versus the Goldstone masses given in Fig. [71 Clearly 
the deviations of the lattice data from unity are consistent with a perturbative expansion. 



Here we use the standard notation, LO = leading order, NLO = next-to-leading order and so on. 
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Quantity 


ATT / ~\ 12 J. 1 

1NLL) nt each process 


NNLO nt 7t^2j^,7P s 


Ci(ir + Y, + ) 


0.66(04)(11) GeV 1 


3.51(18)(25) GeV 1 


C l(7r - ) 


O.o9(0o)(2z) Lrev 


7.44(29) (59) GeV 


Ci(K + v) 


0.44(09) (23) GeV" 1 




C Q1 (if+n) 


0.56(11)(27) GeV" 1 






0.50(06)(14) GeV" 1 




hi 




-0.59(08) (14) GeV" 2 


hl23 




-0.42(10)(10) GeV" 2 



TABLE IV: SU(3) LECs fit from each process at NLO, and from 7T+S+, and tt+S at 
NNLO. The first uncertainty in parentheses is statistical, and the second is the statistical 
and systematic uncertainty added in quadrature. 



Quantity 


LO (fm) 


NLO fit (fm) 


NLO (NNLO fit) (fm) 


NNLO (fm) 




-0.2294 


-0.208(01)(03) 


-0.117(06) (08) 


-0.197(06)(08) 




-0.1158 


-0.105(01) (04) 


0.004(05)(11) 


-0.096(05) (12) 




-0.3971 


-0.311(18) (44) 


0.292(35)(48) 


-0.154(51) (63) 


aRn 


-0.1986 


-0.143(10) (27) 


0.531(28)(68) 


0.128(42)(87) 




-0.4406 


-0.331(12) (31) 


0.324(39)(54) 


-0.127(57)(70) 



TABLE V: SU(3) extrapolated scattering lengths using the LECs from Table IIVI The 
first uncertainty in parentheses is statistical, and the second is the statistical and sys- 
tematic uncertainty added in quadrature. Note that the NLO (NNLO fit) column is using 
d, Cox from the NNLO fit to tt+EV+H . 



At NNLO the situation changes dramatically. This is clear from the plots of T^nlo 
versus the Goldstone masses given in Fig. [71 The shift of the value of V from NLO to NNLO 
is dependent on the renormalization scale //. With the choice /x = A x one would expect this 
shift to be perturbative. However, this is not the case and therefore loop corrections are 
very large at the scale A x . There are many strategies that one may take to fit the LECs in 
the overdetermined system. Here we fit the LECs to the 7r + E + and 7r + S° data, and then 
use these LECs to predict the kaon processes. Therefore, in Fig. [7J, only (a) and (b) are fits. 
The fit LECs are given in Table IIVI While the NNLO LECs h% and appear to be of 
natural size, the NLO LECs Co and Cq% are unnaturally large and therefore are countering 
the large loop effects. The extrapolated 7r + S + and 7r + H° scattering lengths are given in 
Table |V] and appear to be perturbative. Table |V] also gives the extrapolated kaon-baryon 
scattering lengths with the LECs determined from the 7r + S + and 7r + S° data. The resulting 
NNLO predictions deviate by at least 100% from the LO values. Other fitting strategies 
lead to this same conclusion: the kaon-baryon scattering lengths are unstable against chiral 
corrections in the three-flavor chiral expansion, over the range of light-quark masses that we 
consider. 
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FIG. 6: Plots of Tnlo versus the Goldstone masses for the five meson-baryon processes. 
All lattice data is included. 
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FIG. 7: Plots of T NLO and T NNLO versus the Goldstone masses. The line at V — 1 is the 
leading order curve, and dotted line is the physical meson mass. The innermost error bar 
is the statistical uncertainty and the outermost error bar is the statistical and system- 
atic uncertainty added in quadrature. The inner and outer filled bands correspond to the 
statistical and systematic uncertainty, respectively, of the fits to the LECs at NLO and 
NNLO using 7r + S + , and 7r + H° only, for the SU(3) case. 
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VII. SU(2) HBx-PT EXTRAPOLATION 



Given the poor convergence seen in the three-flavor chiral expansion due to the large loop 
corrections, it is natural to consider the two-flavor theory with the strange quark integrated 
out. In this way, 7rS and 7rS may be analyzed in an expansion in m n with no fear of 
corrections that scale as powers of m^. The detailed matching of LECs between the three- 
and two-flavor theories is described in detail in Ref. 18 . We make use of the formulation 



of the 7rE and 7rH T-matrices from 18j| to per form the two-flavor chiral extrapolations for 
O'-k+t.+ i and a^+so. As pointed out in Ref. [18j, there are two representations of the pion- 
hyperon scattering lengths that are equivalent up to omitted higher orders in the chiral 
expansion; one contains a chiral logarithm, and the other is purely a polynomial in m^. 
Using both forms provides a useful check on the systematics of the chiral extrapolation. 



A. Scattering Length Formulas I 



To 0(m%) in the two-flavor chiral expansion, a 7r +s+ and a^+so are given by [1 



1 



47r rn.Tr + mj 



m-. 



2m„ 



2m 2 

H ^r-C 7r + s + + 



J 7T J ', 



rrvt 



m n 2ml 
log 1 — —K 



P 

J 7T 



ft ft n+S ° 2^ ft 



47T + m= 

where the explicit forms — in terms of Lagrangian parameters 
no are given in Ref. 



log h — Ai^+ s o(/x) 

» ft 



(34) 



(35) 



of the LECs C n +^+, h 7r +^+, 
As in the three flavor case, the mixed-action modifi- 
cation to the SU(2) scattering length formula would begin with corrections to the m\ ln(m n ) 



C^+ho and h 
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terms, with the mixed valence-sea pions having the known additive mass shift [66|. We again 
choose fj, = A x = Airfn and evaluate f n at its lattice physical value. In analogy with the 
three-flavor case, we define 



LO 



NLO 



= 1 — C^+BrriT, 



(36) 
(37) 



^nnlo = 1 - Cn+Bmn - h n+B (A x )ml , (38) 

where B is either S + or S°. In Fig. | we give plots of T NLO and T NNLO versus the pion 
mass for the two-flavor case. Clearly the deviations of T from unity are consistent with a 
perturbative expansion at both NLO and NNLO, showing that the loop corrections are much 
smaller at the scale A x than in the three-flavor case. All extracted LECs are of natural size 
and given in Table IVIl The extrapolated 7r + S + and 7r + S° scattering lengths are given in 
Table I VIII The results are consistent with what was found in the three-flavor extrapolation. 
The NLO and NNLO LECs are highly correlated in the NNLO fit. Fig. |9] shows the 68% and 
95% confidence interval error ellipses in the h-C plane for both 7r + S + and 7r + S°. Exploring 
the full 95% confidence interval error ellipse in the h-C plane yields 

a^ +s+ = -0.197 ±0.017 fm ; (39) 
a^ +s0 = -0.098 ± 0.017 fm . (40) 

These are the numbers that we quote as our best determinations of the pion-hyperon scat- 
tering lengths. 
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FIG. 8: T NLO , T NNLO plots for the 7r + £ + , and n + E° processes versus the pion mass. The 
line at T = 1 is the leading order curve, and the dotted line is the physical pion mass. The 
innermost error bar is the statistical uncertainty and the outermost error bar is the statis- 
tical and systematic uncertainty added in quadrature. The inner and outer filled bands 
correspond to the statistical and systematic uncertainty, respectively, of the fits to the 
LECs at NLO and NNLO using 7T+S+, and vr+S for the SU(2) case. 





NLO fit 


NNLO fit 




0.66(04)(11) GeV" 1 


1.98(17)(24) GeV- 1 


Cjt+E 


0.69(06) (22) GeV" 1 


2.01(24) (68) GeV" 1 


K+T.+ 




-0.65(36) (40) GeV" 2 


hn+E 




-0.6(0. 5)(1.1) GeV" 2 



TABLE VI: SU(2) LECs fit from each process at NLO and at NNLO. The first uncer- 
tainty in parentheses is statistical, and the second is the statistical and systematic uncer- 
tainty added in quadrature. 



B. Scattering Length Formulas II 



Ref. [18] makes the interesting observation that replacing f n with its chiral limit value, /, 
yields 



Q-K+T.+ 



27T + ms 



m n .ml ml 

~j2 + + y^ /2 7r+s+ 



Ce + = ^4 + W ;(4i) 



m-- 



47r m n + mn 



m^ ml ml 

~J2 + + -j2 n K+z° 



K+s° = 7^4 + K+s° , (42) 



P 



where t\ is the LEC which governs the pion mass dependence of f v [68|. Note that the 
chiral logs have canceled, and in this form, valid to order ml in the chiral expansion, the 
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Quantity LO (fm) 


NLO (fm) 


NLO (NNLO fit) (fm) 


NNLO (fm) 


a^s -0.2294 
a nE -0.1158 


-0.208(01)(03) 
-0.105(01)(04) 


-0.166(05)(08) 
-0.083(04) (11) 


-0.197(06)(08) 
-0.098(05)(12) 



TABLE VII: SU (2) extrapolated scattering lengths using the LECs from Table I VI I The 
first uncertainty in parentheses is statistical, and the second is the statistical and system- 
atic uncertainty added in quadrature. 



ti + Z + SU(2) 



ji + S° SU(2) 




C v (GeV ) 



FIG. 9: The 68% (light) and 95% (dark) confidence interval error ellipses for fits for the 
7r + S + (left), and tt + E° (right) processes using Eqs. (jMj) and fl35l) . 



scattering lengths have a simple polynomial dependence on m n . Taking the standard value 



/ = 122.9 MeV [18J, |68|] and refitting the LECs yields the results tabulated in Table IVLTT1 
The extrapolated 7r + S + and 7r + S° scattering lengths are given in Table IIXI These results 
are clearly consistent with what was found in the two-flavor extrapolation with the chiral 
logarithm explicit. Fig. [TUl shows the 68% and 95% confidence interval error ellipses in the 
h-C plane for both 7r + S + and 7T + S°. Exploring the full 95% confidence interval error ellipse 
in the h-C plane yields 



0-K+T.+ 
GU+=0 



-0.197 ±0.011 fm ; 
-0.102 ± 0.004 fm . 



(43) 
(44) 



Comparison of these determinations with those of Eq. (140]) give an estimate of the systematic 
error due to truncation of the chiral expansion at order m\. We have also "pruned" the data; 
that is, we have redone all fits omitting the heaviest mass ensemble. While this procedure 
inflates the errors, we see very little shift in the central values. 
In order to plot the scattering length versus m n , we define 



-E+ = 07T+E+ 



1 

2^ 



+ ^C 
+ /2 C 7 



m 
1 



(45) 



(46) 
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NLO fit 



NNLO fit 



CV+S+ 1.28(09)(11) GeV- 1 1.90(10)(17) GeV" 1 

Cy+go 1.84(23)(25) GeV' 1 1.93(12)(48) GeY" 1 
h' x+s+ - -1.33(21)(26) GeV" 2 

feg+go - -1.36(27) (75) GeV' 2 

TABLE VIII: ST/ (2) LECs fit from each process at NLO and at NNLO. The first uncer- 
tainty in parentheses is statistical, and the second is the statistical and systematic uncer- 
tainty added in quadrature. 



Quantity LO (fm) 


NLO (fm) 


NLO (NNLO fit) (fm) 


NNLO (fm) 


a^s -0.2294 

QirE -0.1158 


-0.212(03)(04) 
-0.106(04)(05) 


-0.190(04)(06) 
-0.095(02)(09) 


-0.197(04)(09) 
-0.102(02)(09) 



TABLE IX: SU (2) extrapolated scattering lengths using the LECs from Table IVIIII The 
first uncertainty in parentheses is statistical, and the second is the statistical and system- 
atic uncertainty added in quadrature. 



ti + Z + SU(2) 



7i + E° SU(2) 




C„ v (GeV ) 



C + _o (GeV ) 



FIG. 10: The 68% (light) and 95% (dark) confidence interval error ellipses for fits for the 
7r + S + (left), and (right) processes using Eqs. (j4"Ij) and (|42l) . 



In Fig. [TT] we plot the scattering lengths versus the pion mass. The shaded bands in these 
plots correspond to the standard error in the determination of the LECs, as given in Ta- 
ble Ivnll 

Additional systematic errors arising from the specific lattice formulation that we employ 
are discussed in detail in Ref. and are expected to be well encompassed by our error 
bars. As discussed in section 1111} there is a systematic error in extracting the scattering 
length from the phase shift. We find that range corrections affect the scattering length 
at the 5% level for 7r + S + , and at the 1% level for 7r + S°. Finally, we reiterate that there 
are unquantified systematic errors due to finite- volume and lattice-spacing effects, however, 



23 



these errors are likely encompassed by our quoted errors. 



a xV SU(2) fljt+s o SU(2) 




(MeV) (MeV) 
(a) (b) 

FIG. 11: a plots for the 7r + £ + , and 7r + S° processes versus the pion mass. The diagonal 
line is the leading order curve, and the dotted line is the physical pion mass. The inner- 
most error bar is the statistical uncertainty and the outermost error bar is the statistical 
and systematic uncertainty added in quadrature. The filled bands are the fits to the LECs 
in the SU(2) case at NNLO as in Eqs. (gg) , and (g5J). 



VIII. CONCLUSIONS 

In this paper we have presented the first fully-dynamical lattice QCD calculation of meson- 
baryon scattering. While the phenomenologically most-interesting case of pion-nucleon scat- 
tering involves annihilation diagrams, and therefore, requires more resources than we cur- 
rently have available, we have calculated the ground-state energies of 7r + S + , tt + E°, K + p, 
K + n, and K°E°, which involve no annihilation diagrams. 

An analysis of the scattering lengths of these two-body systems using HB%PT has led 
us to conclude that the three-flavor chiral expansion does not converge over the range of 
light quark masses that we investigate. While the kaon-baryon scattering lengths appear 
perturbative at NLO, a comparison of NNLO with NLO calls into question the convergence 
of the three-flavor chiral expansion. Therefore, we do not quote values for the kaon-baryon 
scattering lengths at the physical point. On the other hand, the 7r + S + and 7r + S° scattering 
lengths appear to have a well-controlled chiral expansion in two-flavor HB%PT. Our results, 
a 7r+£+ = —0.197 ± 0.017 fm, and a n + s o = —0.098 ± 0.017 fm, deviate from the LO (current 
algebra) predictions at the one- and two-sigma level, respectively. We look forward to 
confirmation of these predictions from other lattice QCD calculations and possibly from 
future experiments. 

The HB^PT analyses performed in this work support a general observation about con- 
vergence in the three-flavor chiral expansion, at least for the processes studied here. As the 
pion masses considered in this lattice calculation are comparable to the physical kaon mass, 
the distinct convergence patterns of the two- and three-flavor chiral expansions found in this 
work are suggestive that the breakdown in the three-flavor case is not due to the relative 
largeness of the strange-quark mass as compared to the light quark masses, but rather due 
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to some other enhancement in the coefficients of the loop contributions, possibly related to 
a scaling with powers of n/, the number of flavors. 

While in this paper we have not considered the lowest-lying baryon decuplet, one interest- 
ing process for future study is the ir~Q~ system. It does not involve disconnected diagrams 
since the pions have no valence quarks with the same flavor as the Q~ constituents. It has 
been argued that there is a bound state 69| in this channel, and therefore, it would be of 



interest to determine whether this state appears bound on the lattice at the available quark 
masses. 
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